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ABSTRACT 



o 

o . 

We present the first results of the MIRIAD (MIPS [Multiband Imaging Photometer for Spitzer] 
Infra-Red Imaging of AGB [asymptotic giant branch] Dustshells) project using the Spitzer Space 
Telescope. The primary aim of the project is to probe the material distribution in the extended 
(T) ■ circumstellar envelopes (CSE) of evolved stars and recover the fossil record of their mass loss 

history. Hence, we must map the whole of the CSEs plus the surrounding sky for background 
subtraction, while avoiding the central star that is brighter than the detector saturation limit. 
^ • With our unique mapping strategy, we have achieved better than one MJy sr _1 sensitivity in 

three hours of integration and successfully detected a faint (< 5 MJy sr extended (~ 400") 
far-infrared nebula around the AGB star R Hya. Based on the parabolic structure of the nebula, 
■ the direction of the space motion of the star with respect to the nebula shape, and the presence 

of extended Ha emission co-spatial to the nebula, we suggest that the detected far-IR nebula is 
due to a bow shock at the interface of the interstellar medium and the AGB wind of this moving 
star. This is the first detection of the stellar-wind bow-shock interaction for an AGB star and 
exemplifies the potential of Spitzer as a tool to examine the detailed structure of extended far-IR 
nebulae around bright central sources. 
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1. Introduction 

Evolved stars of intermediate initial mass (~ 1 to 8 M ) are major contributors to the interstellar 
medium (ISM). However, the mechanism by which the ISM is enriched is not well understood. The circum- 
stellar envelopes (CSEs) of evolved stars, asymptotic giant branch (AGB) and post-AGB stars in particular, 
contain the fossil record of their mass loss history, and therefore, have the potential to verify many aspects 
of stellar evolution. IRAS and ISO data indicated that parsec-sized dusty CSEs exist around these infrared 
(IR) objects (e.g. Young, Phillips, & Knapp 1993). Moreover, these large CSEs show evidence for mass-loss 
variations that correlate with evolution of the star (e.g. Speck, Meixner, & Knapp 2000). 

Previous observations lacked the sensitivity and spatial resolution to investigate the full extent and 
structure of these large dusty CSEs. Using the unique sensitivity and mapping capabilities of the Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) on-board the Spitzer Space Telescope (Werner 
et al. 2004), we have been conducting far-IR imaging observations to probe the material distribution in the 
extended CSEs around several evolved stars. Our maps are the deepest yet and have the most complete 
spatial coverage, allowing the most detailed study into the CSE structure to date. 

In this letter, we report the first results of the MIPS IR Imaging of AGB Dustshells (MIRIAD) project 
(Program ID 20258) for the AGB star R Hya. This star is a Mira variable well-known for its decreasing 
period (Zijlstra, Bedding, & Mattei 2002) and the presence of a detached CSE (Young et al. 1993; Hashimoto 
& Izumiura 1997; Hashimoto et al. 1998). Below, we will describe the observations (§2) and data reduction 
(§3), and present the results and discussion (§4). 

2. Observations 

We observed R Hya at 70 and 160^m using Spitzer/MIPS on 2006 February 26 as part of the MIRIAD 
project. R Hya was the first target for which we obtained data in both of the MIPS Ge bands. We mapped 
a roughly square region (24' x 24') at 70/xm and a nearly linear region (24' x 2.'4) at 160/im, while avoiding 
the central region (2' x 1' and 1'.5 x 2' at 70 and 160^im, respectively), using exposures in photometry/fixed- 
clustcr-offset mode. The MIRIAD maps can cover the entire target CSEs plus enough surrounding sky for 
background subtraction. Care was taken to avoid the central star that is brighter than the saturation limit 
of the MIPS arrays. Fig. 1 visualizes the sky coverage. 

The 70/xm observations consisted of three Astronomical Observation Requests (AORs), which sequen- 
tially covered the outer, mid, and inner regions of the CSE. To maximize the dynamic range, we used 
progressively longer exposure times for regions further away from the central star: per-pixel, pcr-Data- 
Collection-Event (DCE) exposure times were 18.87, 62.91, and 251.66 s for the inner, mid, and outer re- 
gions, respectively. The "gap" in the 70/im sky coverage (along the in-scan direction towards the south of 
the central star; Fig. 1) is unavoidable because the central star would fall onto the 24/xm array when the 
70^m aperture covers this region. 

The 160^m offsetting was done by a single AOR so that the linear map was swept in both the forward and 
backward directions along the cross-scan direction. To maximize the dynamic range, we repeated exposures 
at and near both ends of the linear map: per-pixel exposure times are 167.76, 251.56, and 335.52 s depending 
on the number of repeats. Our observing scheme is also summarized in Table 1. 
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3. Data Reduction 

We initially constructed mosaicked maps from the pipeline-calibrated basic calibrated data (BCD) prod- 
ucts using the Mosaicker software (Ver. 093005) provided by the Spitzer Science Center (SSC) 1 . However, 
the BCD-mosaicked maps were affected by detector artifacts. The effects were especially severe in the 70/im 
map, in which the "streaking" (due to residual variations of the slow response of the detector as a function of 
time) was rampant and worsened increasingly towards the end of the observing sequence, probably because 
residual cosmic-ray charge slowly built up in the detector. 

Hence, we performed a custom data reduction from the raw data. First, we used the latest Ge Re- 
processing Tools (GeRT) software (Ver. 041506 of S14 processing) 2 to create our own BCDs. The most 
important step in custom BCD generation is to do proper stimulator-flash (stim) calibration (Gordon et 
al. 2005). For the data from the mid and outer AORs at 70/xm, for which most of the emission is of sky 
(the zodiacal light, ISM, and cosmic IR background), the best results were produced from stim calibrations 
with piece-wise linear fitting excluding stim frames on the nearby point sources. For the data from the 
inner region at 70/im, for which most of the field can be bright by the central star and CSE emission, the 
best outcome was yielded from calibrations with spline fitting using all the raw data. For the 160/xm data, 
for which emission is of both sky and non-sky, the best results were obtained from stim calibrations with 
piece-wise linear fitting using all the raw data. 

For each BCD pixel, we followed the time evolution of pixel values and performed a linear least-squares 
fit using only the "sky" values to determine the baseline for sky subtraction. We tested this baseline fitting 
on the pixel, column, and column-of-four-pixel (readout is done for four columns at a time) bases, and the 
pixel basis fitting yielded the best results. This process worked especially well in removing linear streakings 
in the 70/im maps which resulted from the column-dependent responsivity and the relatively linear way in 
which photometry-mode mapping scans across the sky. 

However, the data obtained from the last offset position at 70/xm showed some residual column- 
dependent artifacts. This probably occurred because the detector responsivity had been compromised by 
the brightest region of the field near the central star scanned immediately before the end of the observing 
sequence (i.e. the detector is due for annealing). The remaining artifacts were removed by flat- fielding all 
the frames taken at the last offset position using their median frame. 

After the baseline fitting, we subtracted the derived sky value at the time of exposure from the entire 
dataset on pixel basis. The subtracted sky values were 13.9, 14.4, and 15.5 MJy sr _1 respectively for the 
outer, mid and inner regions at 70/<m and 19.0 MJy sr" 1 at 160/xm. The increasing sky value with decreasing 
aperture at 70/im illustrates the need of baseline determination far enough away from the central object. 
We adopted 13.9 MJy sr _1 as the 70/im sky emission, and the difference was added back to the BCDs for 
the mid and inner regions (0.5 and 1.6 MJy sr _1 , respectively) to compensate for the over-subtraction. The 
expected sky emission, obtained by the Spot Software 3 , is 11.5 and 14.2 MJy sr" 1 respectively at 70 and 
160/im. This step was done by our own IDL script. 

Finally, the custom-processed BCDs were mosaicked together using the Mosaicker while cosmic-rays 
were removed. The resulting maps still showed some structure intrinsic to the diffraction spikes of the 



1 http: / / ssc.spitzer.caltech.edu/postbcd / 

2 http: / / ssc.spitzer.caltech.edu/mips/gert / 
3 http://ssc. spitzer. caltech.edu/propkit /spot/ 
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point-spread- function (PSF), especially at 70/zm. Thus, we used the Spitzer Tiny Tim (STinyTim) software 4 
to create simulated point-response- function (PRF) maps (with the F v oc v 2 emissivity assumption), and 
subtracted them from the mosaicked maps by scaling flux within an annulus around the central unobserved 
region. 

The resulting maps are in 4'/92 pix -1 and 7'.'99 pix -1 (sub-pixelized from the nominal scale by a factor 
of two) with the average of seven and 15 sky coverages per pixel and the maximum of 45 and 80 sky coverages 
per pixel respectively at 70 and 160^im. At least four sky coverages are achieved in 88% and 94% of the 
mapped area respectively at 70 and 160^.m to fulfill the minimum level of redundancy. The resulting one- 
sigma sensitivities are consistent with the performance estimates done by the SENS-PET tool 5 : we have 
achieved the required sensitivity level of 1 MJy sr _1 as planned. Characteristics of the resulting maps are 
listed in Tabic 2. 



4. Results and Discussion 

Fig. 2 shows the background/PRF-subtracted, mosaicked MIPS color maps of R Hya at 70^m (top 
left) and 160^im (top middle), which clearly display faint, extended emission. Also shown are a pipeline- 
processed ISOPHOT PHT32 map at 60^m (bottom left; Hashimoto & Izumiura 1997) 6 , a HiRes-processed 
IRAS map at lOO^m (bottom middle) 7 , an Ha map (top right) from the Southern H-Alpha Sky Survey 
Atlas 8 (SHASSA; Gaustad et al. 2001) and STinyTim PRF maps at 70 and 160^m (bottom right). 

The 70^m map shows an arc-like surface brightness distribution at > 10 MJy sr _1 surrounded by even 
fainter emission (< 10 MJy sr _1 ) of ~ 400" diameter. The arc structure appears parabolic with the apex 
pointed towards the WNW direction: the distance to the apex from the central star is 100" (see below). While 
the southern part of the arc remains physically thin (~ 40"), the northern part becomes wider (~ 100"). 
The width of the parabola at the tail end is about 370". The brightest region near the apex of the parabolic 
arc (of ~ 20 MJy sr _1 ) can also be seen in the ISOPHOT 60/im map. In retrospect, the ISOPHOT map 
appears to have captured most of the arc to the north, in spite of the instrumental artifacts. 

The 160/im map also displays an extended nebulosity along the arc (< 5 MJy sr _1 ), which appears 
consistent with the 100^m HiRcs/IRAS map as well as the 100/im high-resolution IRAS map based on the 
pyramid maximum entropy method (Hashimoto et al. 1998). There is a bright spot near the NW edge of 
the central unobserved region (at ~ 45" from the star), which does not spatially correspond to the 70/im 
arc. It is probably caused by the near-IR light leak at 160/xm (SSC 2005; For R Hya, mj = —1.3 mag and 
Fi6o/zm/F2^m = 0.0006). 

The integrated flux over the entire nebula is 19.4 ± 3.6 and 5.8 ± 1.3 Jy respectively at 70 and 160^m. 
These values are the lower limits since there is no data in the central unobserved region. From the scaled 
PRFs, the stellar fluxes at 70 and 160^m are estimated to be 33.9 and 4.5 Jy, respectively. 

4 http: / / ssc.spitzer.caltech.edu / archanaly / contributed / 

5 http: / / ssc.spitzcr.caltech.edu / tools / senspet / 

6 http: / / www.iso.vilspa.esa.es/ 

7 http://irsa.ipac. caltech.edu/IRASdocs/hires_over. html 
8 http://amundscn. swart hmore.edu/SHASSA/ 
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Van Buren & McCray (1988) have found parabolic far-IR nebulae around hot stars and attributed some 
of them to bow shocks at the interface between the stellar wind and ISM around a moving star. R Hya is 
known to have the proper motion of (-60.73, 11.01) mas yr -1 (Perryman et al. 1997), towards the apex of 
the arc as indicated by the bottom left arrow in Fig. 2. Moreover, the shape of the far-IR arc closely follows 
the curve y = x 2 /3l (where / is the distance between the apex and star), which has been shown numerically 
to represent the apex shape of the stellar- wind bow-shock by Mac Low et al. (1991), as delineated by the 
dashed lines in Fig. 2. Furthermore, Ha emission is extended along the parabolic apex and axis, showing 
the presence of ionized gas in these regions. 

Therefore, we attribute the parabolic far-IR nebula around R Hya to shock-excited line emission (such as 
[O I] 63 & 146^m) and locally heated dust emission arising from a stellar-wind bow-shock interface between 
the ISM and the swept-up AGB wind of this moving star. This is the first detection of the stellar-wind 
bow-shock interaction for an AGB star, and is a surprise discovery given our original aim of the study. 
Nevertheless, our results exemplify the potential of Spitzer as a tool to examine the far-IR structure of 
extended emission around the bright central object. A follow-up IR spectroscopy can determine the emission 
characteristics. 

At 165 pc (Zijlstra et al. 2002) with the -10 km s _1 radial velocity (Knapp et al. 1998), the space 
velocity of R Hya is 50 ± 1 km s _1 into the direction 12° away (receding side) from the plane of the sky. 
The asymmetric CO profiles (Knapp et al. 1998; Teyssier et al. 2006) can then be explained by this inclined 
bow-shock, in which the receding and approaching shock fronts intersect with the line of sight at distinct 
angles. 

The physical dimensions of the observed bow shock (0.03 pc thick located at 7.9 x 10~ 2 pc ahead of the 
star) are consistent with the stellar- wind bow-shock models for AGB stars (Wareing et al. 2006a) . Adopting 
the mass-loss rate of 3 x 10~ 7 M Q yr -1 (Zijlstra et al. 2002), wind velocity of 10 km s _1 (e.g., Knapp et al. 
1998), and formula for I derived from the momentum conservation across the shock (eq.2 of Van Buren et 
al. 1990), the ambient H density is 0.4 cm~ 3 and the shock surface density is 2.5 x 10~ 7 g cm~ 2 . Thus, the 
estimated amount of matter contained in the bow-shock is about 1.3 x 10~ 4 M Q , assuming a paraboloidal 
bow-shock of y = x 2 /3l with a 300" opening radius. 

The MIPS maps indicate the presence of matter in the downstream (emission at ~ 3 MJy sr _1 up to 
300" away, which is corroborated by the Ha map). The downstream emission may arise from a tail of a 
steady-state bow shock that is formed by ram-pressure stripping from the head of the shock (Wareing et 
al. 2006a). The local emission peaks at <~ 250" away in the downstream may correspond to vortices that 
results from instabilities in the downstream flow shed off from the bow shock (C. J. Wareing, in preparation). 
Modeling of the bow shock nebula of R Hya will be presented by Wareing et al. (2006b). 

This work is based on observations made with the Spitzer Space Telescope, which is operated by the 
JPL/Caltech under a contract with NASA. Support for this work was provided by NASA through an award 
issued by JPL/Caltech. We also acknowledge additional support for the following individuals: for Ueta 
by an NPP Research Fellowship Award, for Speck by NASA ADP grant (NAG 5-12675), for Herwig by the 
LDRD program (20060357ER) at LANL, for Gchrz by NASA (Contract 1215746) issued by JPL/Caltech, for 
Szczerba by grant 2.P03D. 017.25, for Izumiura by Grant-in-Aid (C) from JSPS (No.17540221), for Matsuura 
by JSPS, and for Elitzur by NSF AST grant (0507421). We thank C. J. Wareing for sharing his insights on 
the stellar-wind bow-shocks with us. 
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Table 1. Summary of Spitzer/MIRIAD Observations of R Hya 



Image Scale Field Exposure Integration 3, DCE Total Duration 

Band AOR KEY Start Time (arcscc) Size (sec) Cycle (sec) Offsets Counts (sec) 



70nm 


(outer) 


14472448 


2006 


Feb 


26 


07:15:03 


9.84 


Large 


10 


2 


125.83 


13 


418 


6493 


70nm 


(mid) 


14472704 


2006 


Feb 


26 


08:59:52 


9.84 


Large 


10 


1 


62.91 


8 


145 


2625 




(inner) 


14472960 


2006 


Feb 


26 


09:40:12 


9.84 


Large 


3 


1 


18.87 


6 


109 


1204 


160tir 




14472192 


2006 


Feb 


26 


04:20:36 


16.0 


Small 


10 


2 


83.88 


22 


792 


10256 



integration times are per-pixcl pcr-DCE in "real" second. 
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Tabic 2. Summary of Measurements and Map Characteristics 



Band 


Flux (Jy) 






Sky, a 


(MJy Br" 1 ) 




Coyer. 


rgc (pix- 1 ) 


CSE a 


Star D 


Sky 




"best 


CT cst c 


Avg 


Best 


70/xm (mosaic) 


19-4 ± 3.6 


33.9 


13.9 


0.47 


0.15 




7 


45 


70^m (outer) 






13.9 


0.56 


0.16 


0.26 


5 


45 


70(»m (mid) 






14.4 


0.68 


0.21 


0.37 


3 


24 


70/im (inner) 






15.5 


0.86 


0.28 


0.68 


2 


17 


160>im 


5.8 ± 1.3 


4.5 


19.0 


0.31 


0.12 


0.38 


15 


80 



a Lowcr-limit. 

b Bascd on the scaled PSF/PRF maps. 

c One-sigma sensitivity estimated with the SENS-PET tool. 
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Fig. 1. — MIRIAD sky coverage over R Hya: at 160/im (thick white line) and at 70/Ltm for the outer (thin 
white line), intermediate (black line), and inner (gray line) regions. The "gap" to the south of the star is 
unavoidable to prevent the central star from falling into the 24/itm array. The central box (dashed white 
line) delineates the field of view of Fig. 2. The grayscale SERCT plate image (at 790 nm) was taken from 
the Digitized Sky Surveys. 
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Fig. 2. — The background/PRF-subtracted, mosaicked MIPS maps of R Hya (N is up, E to the left) at 70/im 
(top left) and 160/xm (top middle). For comparison, an ISOPHOT PHT32 map at 60/im (bottom left), a 
HiRes-processed IRAS map at lOO^m (bottom middle), a SHASSA Ha map (top right), and STinyTim 
PRF maps (bottom right) are also shown. The images are zeroed at the position of R Hya ( [a, <5]2ooo — 
[13 h 29 m 42!7803, -23°16'52"792]; Perryman et al. 1997). The position of the star is indicated by the "star". 
The tickmarks indicate the angular offsets in arcsec. Linear color-scaling of surface brightness (in MJy sr _1 ) 
is shown at the bottom of each panel. A parabolic curve y = x 2 /3l, which closely represents the stellar-wind 
bow-shock near the apex, is displayed by the dashed lines. The arrow at the bottom left corner shows the 
direction of the proper motion of the star, (-60.73, 11.01) mas yr~ x . 



